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IGR J16351-5806: another close by Compton-thick AGN 
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ABSTRACT 



1 INTRODUCTION 



Compton thick AGN are by definition those in which the X-ray 
obscuring matter has a column density equal to, or larger than, 
the inverse of the Thomson cross-section (N w > cr^. 1 =1.5 x 10 24 
cm" 2 ). The cross-sections for Compton scattering and photoelec- 
tric absorption have approximately the same value for energies of 
the order of 10 keV, which incidentally coincides with the lower 
energy threshold of hard X-ray telescopes like INTEGRAL/IBIS 
and Swift/BAT. For columns in the range 10 24 -10 25 cm" 2 the nu- 
clear radiation is still visible above 10 keV and the source is called 
mildly Compton thick. For higher values of N H , the entire high 
energy spectrum is down-scattered by Compton recoil and hence 
depressed over the entire X/gamma-ray band. Generally indirect 
arguments allow determination of the Compton thick nature of an 
AGN, such as the presence of a strong iron K line complex at 6.4 
- 7 keV and the characteristic reflection spectrum. The presence 
of Compton thick matter can also be inferred through the ratio of 



IGR J16351-5806 has been associated with the Seyfert 2 galaxy ESO 137-G34, hav- 
ing been first reported as a high energy emitter in the third INTEGRAL/IBIS survey. Us- 
ing a new diagnostic tool based on X-ray column density measurements vs softness ratios 
(F2-10 yfc e y/F2o-ioo kev), Malizia et al. (2007) identified this source as a candidate Compton 
thick AGN. In the present work we have analysed combined XMM-Newton and INTEGRAL 
data of IGR J16351-5806 in order to study its broad band spectrum and investigate its Comp- 
ton thick nature. The prominent Ka fluorescence line around 6.4 keV (EW > 1 keV) together 
with a flat 2-10 keV spectrum immediately point to a highly obscured source. The overall 
spectrum can be interpreted in terms of a transmission scenario where some of the high en- 
ergy radiation is able to penetrate through the thick absorption and be observed together with 
its reflection from the surface of the torus. A good fit is also obtained using a pure reflection 
spectrum; in this case the primary continuum is totally depressed and only its reflection is 
observed. An alternative possibility is that of a complex absorption, where two layers of ab- 
sorbing matter each partially covering the central nucleus are present in IGR J16351-5806. 
All three scenarios are compatible from a statistical viewpoint and provide reasonable AGN 
spectral parameters; more importantly all point to a source with an absorbing column greater 
than 1.5 x 10 24 cm" 2 , i.e. to a Compton thick AGN. Because of this heavy obscuration, some 
extra components which would otherwise be hidden are able to emerge at low energies and 
can be studied; a thermal component with kT in the range 0.6-0.7 keV and free metal abun- 
dance is statistically requied in all three scenarios while a scattered power law is only present 
in the pure reflection model. By providing strong evidence for the Compton thick nature of 
IGR J16351-5806, we indirectly confirm the validity of the Malizia et al. diagnostic diagram. 
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isotropic versus anisotropic emission; primarily the Lx/L[o///j raticQ 
has been used in the literature to recognise Compton thick sources 
(Bassani et al. 1999 and Panessa and Bassani 2002). There are var- 
ious reasons why these AGN are important observational targets. 
They are a key ingredient of the Cosmic X-ray Background (CXB) 
and, through their absorbed/re-emitted radiation, also of the IR one. 
The heavy obscuration present in these objects allows some spec- 
tral components which are otherwise usually hidden by the nuclear 
emission to emerge and be studied in some detail. Despite their 
importance, studies of Compton thick AGN have so far been ham- 
pered by the lack of high energy coverage in X-ray missions and 
subsequent difficulties in performing surveys above 10 keV. This 
situation has been changed by the advent of satellites like INTE- 
GRAL and Swift, which are surveying a large fraction of the sky 
with good sensitivity and arcminute location. As a result, a num- 



1 Where the L[o;;/j is corrected for reddening in the host galaxy by means 
of the H a /Hp B aimer decrement 
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ber of Compton thick candidates are now emerging and can be 
the subject of more in-depth studies (Ueda et al. 2007, Comas- 
tri et al. 2007, Sazonov et al. 2008). In particular, Malizia et al. 
(2007) studying a sample of 34 INTEGRAL AGN in the X-ray 
band, proposed a new diagnostic tool which allows the isolation 
of Compton thick AGN candidates: this relies on measurement 
of the source X-ray column density and softness ratio defined as 
F2-10 kev[FiQ-\oo kev- Using this tool, Malizia and co-workers sug- 
gested 3 INTEGRAL-detected AGN as Compton thick candidates 
(IGR J14175-4641, IGR J16351-5806 and Swift J0601.9-8636) 
based on the apparent absence of an absorbing column density ac- 
companied by a very low softness ratio. One of these sources, Swift 
J0601. 9-8636, has subsequently been observed by Suzaku (Ueda 
et al. 2007) and found to be a Compton thick AGN: the Suzaku 
spectrum revealed a heavily absorbed power law component with a 
column density of N H ~ 10 24 cirr 2 and an intense reflection com- 
ponent with a solid angle > 27r from a cold optically thick medium. 
In this paper we present XMM-Newton and INTEGRAL data of IGR 
J 1635 1-5806 which provide strong evidence for the Compton-thick 
nature of also this object and indirectly confirm the validity of the 
diagnostic diagram proposed by Malizia et al. 



2 IGR J16351-5806 = ESO 137-G34 

IGR J 1635 1-5806 was first reported as a hard X-ray emitter by Bird 
et al. (2007) and subsequently confirmed by Krivonos et al. (2007). 
Here and in the following, we use data collected for the 3rd IBIS 
survey (Bird et al. 2007), which in this region correspond to a total 
exposure of ~ 1.4 Ms. A clear excess is present in the 18-60 keV 
IBIS map with a significance of ~8.3cr at a position correspond- 
ing to RA(2000)=16 h 35 m 11.04 s and Dec(2000)=-58°05'24.0" and 
with an associated uncertainty of 3.1' (90% confidence level). Our 
position is compatible with that reported by Krivonos et al (2007) 
(note that their 2.1' error radius corresponds to 68% confidence 
level). The source position was then reduced to arcsec accuracy 
by means of Swi/f/XRT observations which located the source 
atRA(J2000)=16 h 35 m 13.17 s ,Dec(J2000)=-58°04'49.68" (5.1 arc- 
sec uncertainty) and confirmed the association with the Seyfert 2 
galaxy ESO 137-34 at z = 0.009 (Malizia et al. 2007, Landi et al. 
2007). 

ESO 137-G34 is an SO/a galaxy which in HST (Hubble Space Tele- 
scope) images displays a long dust lane east of the nucleus, running 
northwest-southeast, parallel to the photometric major axis of the 
galaxy and a second dust lane perpendicular to the first one (Fer- 
ruit et al. 2000); moreover in the HST [O III]/[N II] + Ha map, the 
morphology of the central regions appears biconical. 
The source was observed in the radio waveband with the ATCA 
telescope at 3.5cm: it has a flux of 9.6 mjy and displays a com- 
plex morphology consisting of three knots aligned along a position 
angle of -40 degrees (Morganti et al. 1999). The radio emission is 
coincident with the inner part of the line-emitting gas. 
The source has been studied for the first time in X-rays by Malizia 
et al. (2007) using Swift/XRT observations, and its X-ray spectrum 
was broadly defined as being a flat (r= 1.6) power law with no ab- 
sorption in excess of the Galactic value. The 2-10 keV observed 
flux is 3 x 10~ 13 erg cirr 2 s _I almost a factor of 100 weaker than 
the 20-100 keV flux measured by INTEGRAL (2 x 10"" erg crrT 2 
s~'). IGR J116351-5806 was serendipitously observed by XMM- 
Newton on February 13th, 2006 and these X-ray data together with 
the INTEGRAL ones are used in this work to study the broad band 




Figure 1. XMM-Newton pn 4-10 keV image of IGR J 1635 1-5 806 with the 
3.1 arcmin INTEGRAL error box superimposed 

(0.5-100 keV) spectrum of the source. This analysis provides evi- 
dence for the presence of features typical of a Compton thick AGN. 



3 DATA ANALYSIS 

XMM-Newton EPIC-pn (Struder et al. 2001) data were reprocessed 
using the XMM Standard Analysis Software (SAS) version 7.1.2 
employing the latest available calibration files. Only patterns cor- 
responding to single and double events (PATTERN<4) were taken 
into account; the standard selection filter FLAG=0 was applied. 
The observation has been filtered for periods of high background 
and the resulting exposure is 17 ks. Source counts were extracted 
from circular region of radius 25" centred on the source in order 
to exclude the extended emission associated to the galaxy; back- 
ground spectra were extracted from circular regions close to the 
source, or from source-free regions of 80" radius. The ancillary re- 
sponse matrices ( ARFs) and the detector response matrices (RMFs) 
were generated using the XMM SAS tasks arfgen and rmfgen; spec- 
tral channels were rebinned in order to achieve a minimum of 20 
counts for each bin. Figure 1 shows the 4-10 keV pn image of the 
region surrounding IGR J1635 1-5806 with the INTEGRAL error 
box superimposed; it clearly shows that the AGN is the only source 
of hard X-ray photons and hence the counterpart of the INTEGRAL 
source. Inspection of the XMM light curve provides no evidence 
for variability over the observing period, so all data were used for 
the spectral analysis. 

The INTEGRAL data reported here consist of several pointings per- 
formed by the IBIS/ISGRI (Ubertini et al. 2003, Lebrun et al. 2003) 
instrument between revolution 12 and 429, i.e. the period from 
launch to the end of April 2006. ISGRI images for each available 
pointing were generated in various energy bands using the ISDC 
offline scientific analysis software OSA (Goldwurm et al. 2003) 
version 5.1. Count rates at the position of the source were extracted 
from individual images in order to provide light curves in various 
energy bands; from these light curves, average fluxes were then es- 
timated and combined to produce an average source spectrum in 
the 20-110 keV band (see Bird et al. 2007 for details). 
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Figure 2. X-ray and Gamma ray data compared to the simple power law 
used to fit the 2-10 keV data (see text to details); low energy data are from 
XMM-Newton EPIC-pn instrument while high energy points are from IN- 
TEGRALIBIS/ISGRI telescope 

The XMM-Newton data have been fitted together with the IN- 
TEGRAL ones using XSPEC v. 12.3. 1 ; errors are quoted at 90% con- 
fidence level for one parameter of interest (A^ 2 =2.71). In the fit- 
ting procedure, a cross-calibration constant, C, between X-ray and 
gamma-ray data has been introduced and fixed to 1 since gener- 
ally a good agreement between the two instruments has been found 
in previous works regarding various AGNs (De Rosa et al. 2008, 
Panessa et al. 2008, Molina et al. 2008). In the following analysis, 
we always include the Galactic column density which in the direc- 
tion of IGR J1635 1-5806 is N# Goi =2.47 x 10 21 cm~ 2 (Dickey & 
Lockman 1990). 

A simple absorbed power law model fails to reproduce the broad 
band spectrum of IGR J 1635 1-5806: it gives a steep photon index 
(r=2.5), no absorption and a x~ °f 607. 1 1 for 93 d.o.f. However 
if only the XMM data are considered over the 2-10 keV band, the 
photon index is much flatter (r <sc 1), there still is no evidence of 
absorption in excess to the Galactic value, and the observed 2-10 
keV flux is ~ 5 x 10~ 13 erg cnT 2 s -1 . This spectrum is similar to 
that first observed in the prototype Compton thick AGN NGC 1068 
(Matt et al. 2000). Although the XMM spectrum is flatter than the 
XRT one (probably due to the different energy band used), the X- 
ray fluxes are compatible, suggesting no significant variation over 
one year time scale. The use of the XMM fit over a broader band 
(0.2-100 keV, see figure 2) helps in identifying extra features which 
contribute to the overall emission. Indeed, inspection of figure 2 in- 
dicates the presence of a strong soft excess and a prominent line at 
around 6.4 keV. These features, together with the flat photon index, 
strongly resemble a Compton thick AGN. 

The shape of the emission below 2 keV and the presence of emis- 
sion lines clearly indicate the presence of gas photoionized by 
the central engine of the AGN (see discussion). This component 
is often present in Compton thick AGN (Guainazzi et al. 1999, 
Guainazzi et al. 2005) where, thanks to the heavy absorption, it 
emerges and can be well studied. It is generally fitted in XSPEC with 
amekal model having free metal abundance. Indeed the addition of 
a mekal model improves the fit (Ajf 2 =377 for 2 d.o.f.) and gives a 
gas temperature kT of 0.64*2'?; keV and, as already found in other 
absorbed AGN, a low metal abundance value A z of 0.07*2™, still 
leaving an extremely flat power law continuum (r~0.3). Although 
this model yields an adequate fit to our spectrum, it is worth not- 
ing that it may be an oversimplified parameterization of the data. 
High resolution spectroscopy of nearby Seyfert 2s have, in fact, 
demonstrated that the soft X-ray emission is often dominated by 
emission lines with negligible contribution by an underlying con- 



tinuum. Blending of these emission lines in the EPIC spectra can 
mimic a continuum emission (Iwasawa et al. 2003). As our main 
goal is to prove the Compton thick nature of IGR J1635 1-5806, the 
uncertainties induced by a purely phenomenological modelling of 
the soft excess will not substantially affect the core results of our 
paper. 

Another significant improvement (Ajf 2 = 103 for 3 d.o.f.) is obtained 
when we introduce the Ka iron fluorescence emission line: the en- 
ergy of the line is found to be at 6.39 ± 0.02 keV with an equiv- 
alent width of about 2 keV. The intrinsic width of the line is nar- 
row (0.09±0.04) and for simplicity has been frozen to the observed 
value in subsequent fits. The flat slope of the power law (~0.3) as 
well as the strength of the line at 6.4 keV (EW > 1 keV) and the 
lack of absorption again suggest that the source is highly absorbed 
and likely to be in the Compton thick regime (Matt et al. 2000). 
Two models are generally used to account for the soft gamma-ray 
continuum of AGN in this regime: 

(i) a transmission model more appropriate to mildly Compton 
thick AGN; in this scenario the power law is absorbed by a column 
density in the range 10 24 - < 10 25 and a Compton reflection com- 
ponent is often present and added to the fit (model MEKAL + WA * 
PO + ZGA + PEXRAV in XSPEC0 

(ii) a pure reflection model more typical of Compton thick AGN 
with N// > 10 25 cirr 2 ; in this case the absorbed power law is no 
longer detected as it is totally absorbed and only the reflection com- 
ponent is visible (model MEKAL + PEXRAV + ZGA in XSPEC). 

In order to evaluate the contribution of the Compton reflection 
component, we used the code developed by Magdziarz & Zdziarski 
(1995), where we fixed the cosine of the inclination angle to 0.5 and 
the cut-off energy at 1 MeV. 

In the transmission scenario we find a column density of N H = 
4.8 xlO 24 cirr 2 , a photon index of 1.8 and a reflection of ~1 
0f 2 =98.6/86); it is worth noting that in this scenario the reflection 
component is highly requested as without it the model is unable to 
reproduce the INTEGRAL data. 

The reflection model yields a good fit too (89.45/88): the high en- 
ergy part of the spectrum is now represented by a reflection of 0.6 
while the column density is indirectely estimated to be N# > 10 25 
cnr 2 , i.e. a value typical of reflection dominated Seyfert 2 galax- 
ies. In both models, the temperature and abundance of the thermal 
component are consistent while the iron line mantains an equivalent 
width around 1.5 keV. Both scenarios are able to explain such high 
values of EW and are therefore equally viable (Ghisellini, Haardt 
& Matt 1994). Despite the good fits, neither model (transmission 
or reflection) accounts for the excess counts in the range 2-5 keV 
which are still visible in the residuals and require the introduction 
of an additional component in the form of a power law. This ad- 
dition is marginally required by the data in the transmission model 
(~90%), while it is more significant in the reflection model (>99%). 
If we assume that this extra component is due to scattered emission 
from the primary continuum, then the scattered fraction is around 
6% for the reflection model and 10% for the transmission one; both 
values are similar to those typically observed in Seyfert 2 galaxies. 
If present, this soft component may be associated with the radio 
knots and emission line cones observed by ATCA and HST, or it 
could be due to soft emission from the galaxy and/or starburst ac- 
tivity in the source. 

2 this is the model also adopted by Ueda et al. (2007) to describe the Suzaku 
spectrum of Swift J0601.9-8636 
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Figure 3. X-ray (XMM-Newton) and y-ray (INTEGRAL-IBIS) data fitted 
with the complex absorption model (mekal+pcfabs*pcfabs*po+zga) 



Note however that the addition of this extra component in the trans- 
mission model flattens the primary continuum (T = 1.31^JHS), gives 
a smaller column density (N w ~ 9x 10 23 cm -2 ) and has more dif- 
ficulties in explaining a large iron line EW. On the basis of these 
indications, and in view of the marginal evidence in the data of 
this extra component, we prefer to neglect it in the transmission 
scenario but keep it in the reflection one; the relative best fit param- 
eters for both models are given in columns 1 and 2 of table 1 . 
We also consider an alternative scenario where IGR J 1635 1-5806 is 
characterised by a complex absorption partially covering the central 
nucleus ( MEKAL + PCFABS * PCFABS * PO + ZGA in XSPEC). 
This model mimics the reflection component and provides a canon- 
ical AGN spectrum (r=1.8) for the source (see table 1, column 3). 
The absorption required by the data is in the form of two columns 
(N H] ~ 4 x 10 24 citT 2 and N„ 2 ~ 6 x 10 23 citT 2 ), both covering 80- 
90% of the source. Also in this case the combinations of such high 
columns are able to explain the observed iron line EW (Ghisellini, 
Haardt & Matt 1994). 

Despite being rarely used in the literature to model Compton thick 
AGNs, this model is of interest in view of recent studies on the 
torus geometry and nature which strongly indicate that this struc- 
ture is clumpy and made of dusty clouds that are individually opti- 
cally thick (Elitzur 2008). It is possible that in IGR J16351-5806, 
the torus is stratified with a Compton thick region which progres- 
sively decreases in column density moving away from its plane; in 
this case two layers with different absorption (one Compton thick 
and one Compton thin) are intercepted by our line of sight (Risaliti 
et al. 2005); if this is the case, variations in the column densities 
are expected as the clouds move through these two regions. Be- 
sides giving the best x 1 (88.7/85) of all 3 models used here, the 
complex absorption scenario is also attractive for its simplicity as 
it does not require additional spectral features such as the reflec- 
tion and the extra power law components. Although further X-ray 
observations of the source can confirm or deny the presence of a 
complex absorber by measuring variations in the column densities, 
we tentatively assume this to be the more appropriate scenario for 
this source. The unfolded spectrum of IGR 116351-5806 fitted with 
this model is shown in figure 3. 

4 DISCUSSION AND CONCLUSIONS 

In this work we have analysed the combined XMM-Newton and 
INTEGRAL-IBIS data of the Seyfert 2 galaxy IGR J16351-5806. 
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transmission reflection^' complex abs f 



kT 


n 69+0.06 
-Q04 


66 +0 06 


6? +0 05 

u '° -0.05 


A z 




031^ 


>0.06 


N«, (cm- 2 ) 


4.8+™ x 10 


24 


3.7 +4 -* x 10 24 








N„ 2 (cm- 2 ) 






5.8 +22 x 10 23 


Cj5 






n q 5+0.08 
• -0.09 


*soft 




242+0.42 




R 

E lme (keV) 


1 0-1+0.21 
-0.82 

6 39+ 8 - 02 


i tt+o:89 

-0.07 
0,59+0-83 

6 39+8S 


1 7O+0.08 

6 39+° 02 

o - Jy -0.02 


EW (keV) 


1 5 +0 - 27 
i.j_ 022 


1 56+ - 28 
1 - J -0.29 


1 36+ 050 

-0.38 


* 2 /dof 


98.59/86 


89.45/86 


88.70/85 


F2-10faV 




~ 5 x 1(T 13 cm" 2 s" 


1 


F20-100ieV 




~2x 10-' 1 cm" 2 s" 


l 



(fl) MEKAL + WA * P0 + ZGA + PEXRAV 

(W MEKAL + P0 + PEXRAV + ZGA 

(c) MEKAL + PCFABS * PCFABS * PO + ZGA 



The broad band spectrum of IGR J 1635 1-5806 provides immedi- 
ate strong indications of its Compton thick nature, in particular the 
presence of a flat 2-10 keV spectrum and of a prominent emission 
line around 6.4 keV. The overall spectrum can either be interpreted 
in terms of a transmission scenario where some of the high energy 
radiation is able to penetrate through the thick absorption and be 
observed together with its reflection from the surface of the torus. 
A good fit is also obtained using a pure reflection spectrum: in this 
case the primary continuum is totally depressed and only its reflec- 
tion is observed. An alternative possibility is that of complex ab- 
sorption, where two layers of absorbing matter each partially cov- 
ering the central nucleus, are present in IGR J1635 1-5806. All three 
scenarios are compatible with the data from a statistical view point, 
provide reasonable AGN spectral parameters, and are capable of 
producing the observed iron line EW. They are therefore indistin- 
guishable using the present data set, although the complex absorp- 
tion model is more attractive for its simplicity and its geometrical 
implications for the torus that are more in line with recent obser- 
vational and theoretical findings. Whether mildly absorbed, totally 
hidden or partially covered, the source has a column density ex- 
ceeding 1.5 xlO 24 cmr 2 , which confirms IGR J1635 1-5806 as a 
new Compton thick AGN detected by INTEGRAL in our neigh- 
bourhood. 

Because of this heavy obscuration, some extra components which 
would be otherwise hidden, are able to emerge at low energies and 
can be studied; this includes the presence of thermal emission with 
free metal abundance and possibly of a scattered power law. The 
presence in the soft X-ray band of a hot diffuse gas with kT ~0.6- 
0.7 and A z >0.05 is in agreement with observational evidence com- 
ing from grating instruments onboard Chandra and XMM-Newton 
on other Compton thick Seyfert 2 galaxies (Guainazzi et a. 2005): 
within the sample analysed by these authors similar temperature 
and abundance values were found and interpreted in terms of emis- 
sion from an optically thin, collisionally ionized plasma. Although 
this could also be the interpretation in the case of IGR J 16351- 
5806, we note the model used to fit the soft part of the spectrum of 
IGR 116351-5806 is still oversimplified and we need RGS data for 
a more detailed analysis; unfortunately the source is too weak for 
good quality RGS data in the present observation and needs to be 
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re-observed with a longer exposure. 

IGR J16351-5806 has been highlighted previously by Malizia et 
al. (2007) as one of three possible Compton thick objects, using a 
new diagnostic tool based on the measurement of the X-ray column 
density and softness ratio F2-10 fa-v/F20-ioo tev- Of their three candi- 
dates, SWIFT J0601. 9-8636, has been subsequently confirmed as 
a Compton thick AGN thanks to Suzaku observations (Ueda et al. 
2007), another, IGR J14175-4641 is very dim and with too lim- 
ited spectral information below 10 keV to be properly studied - 
we hope to overcome this with a future XMM pointing. The third 
source is IGR J16351-5806 discussed here. By proving its Comp- 
ton thick nature we have also indirectly demonstrated the validity 
of the Malizia et al. (2007) diagnostic diagram and provided an 
extra tool in the search for such rare AGN. The use of this diag- 
nostic together with other information will be further tested on a 
larger INTEGRAL AGN sample for which 2-10 keV data are now 
becoming available. 
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